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We report on an updated Paris NN optical potential. The long- and intermediate-range real parts 
are obtained by G-parity transformation of the Paris NN potential based on a theoretical dispersion- 
relation treatment of the correlated and uncorrelated two-pion exchange. The short-range imaginary 
potential parametrization results from the calculation of the NN annihilation box diagram into two 
mesons with a nucleon-antinucleon intermediate state in the crossed channel. The parametrized real 
and imaginary short range parts are determined by fitting not only the existing experimental data 
included in the 1999 version of the Paris NN potential, but also the recent antiprotonic-hydrogen 
data and np total cross sections. The description of these new observables is improved. Only this 
readjusted potential generates an isospin zero ^So, 52 MeV broad quasibound state at 4.8 MeV 
below the threshold. Recent BES data on J/tp decays could support the existence of such a state. 
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1. INTRODUCTION 

There has been recently a renewal of interest in the 
nucleon-antinucleon, NN, interaction due to the obser- 
vation of near threshold enhancements in the proton- 
antiproton, pp, invariant mass spectrum of heavy meson 
decays such as J/ip ^pp B'^ ppK'^, 
ppn+ 0,5°^ ppK*° and B° D^*^^pP |4i]. On the 
other hand, no such structure was observed by the BES 
Collaboration for the J/ip ^ n^pp decays For the 
radiative and pionic J/ tp decays reported in Ref . [l| , two 
of us have proposed a natural explanation following from 
a traditional model of pp interactions Q. These inter- 
actions originate from the Paris NN potentials 
Taking into account the low energy allowed final states, 
the BES data are well reproduced with an isospin one 
ppi^Pi) wave for the ir'^pp channel and a pp(^S'o) wave 
for the ^pp channel. It was furthermore shown in Ref. [5] 
that the best results were obtained with an upgraded 
NN Paris potential constrained not only by the set of 
data used in the 1999 version [sj, but also by recent to- 
tal np cross sections of Ref. [9] and antiproton-hydrogen 
widths and shifts [l3j[IH- Only this recently readjusted 
potential has an isospin T = Sq quasibound state close 
to the pp threshold. The existence of such a state has 
some support from the BES data even if the low-energy 
pp spectrum of the radiative decay could also be repro- 
duced in Ref. [l^l using the T = 1 S'-wave of the meson- 
exchange Jiilich-Bonn NN model where no "'^S'o bound 
state is present. The aim of the present work is to report 
on the updated Paris NN potential used in Ref. [5]. 

The paper is organized as follows. In Sec. [H after a 
brief reminder of the model, we compare its results to 
the experimental scattering observables and to the re- 
sults of the 1999 NN Paris potential. This comparison is 
also done for the antiprotonic-hydrogen level shifts and 



TABLE I: Heights of the different real potentials U(r) at 
r = ra = 0.188 fm and r = r2 = 0.587 fm together with 
the parameters g (dimensionless) and / of the imaginary po- 
tentials. These quantities, determined by the fit to experi- 
mental observables, are compared with those of the Paris 99 
potential All U{r) are in MeV but the U^{r) which are 
dimensionless. The definitions of the real and imaginary po- 
tentials can be found in the Appendix. 
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widths. The close to threshold bound states and reso- 
nances of the present work are searched for and compared 
to those of the 1999 model. We then plot and compare 
the optical potentials of both solutions. Section [3] is de- 
voted to a discussion of our results. A summary and con- 
clusions are presented in Sec. HI Finally, the Appendix 
reminds the reader of the full expression of the optical 
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FIG. 1: Total and annihilation cross sections for the pp and np systems. The references of the experimental data can be found 
in Ref. The data of lazzi et al. in (c) are from Ref. jSli]. 



potential. 

2. MODEL AND RESULTS 

2.1. Brief reminder of the model 

The 1982 Paris NN opticalpotential @ was itself read- 
justed in 1994 and 1999_[8'l. For all these potentials 
and the present one, the NN interaction is described by 
an energy dependent optical potential 

Vj^N (r, TLab) = Unn (r, TLab) - i W^n (r, TLab) , (1) 

where the nonlocality of the real C/^^ and imaginary 
Wj^fj^ potentials are accounted for by a linear energy de- 
pendence in the kinetic energy T^ab- Meson exchanges 



explain in a satisfactory way the NN force for large and 
medium distances between the nucleons. Therefore the 
long and intermediate range real parts, i.e. those for 
inter NN distances r > 1 fm, are obtained by the G- 
parity transformation of the corresponding parts of the 
Paris NN potential [l^ based on a theoretical dispersion- 
relation treatment of the correlated and uncorrelated 
two-pion exchange These real potentials contain, 

besides the one-pion exchange, the two-pion exchange 
and the oj and Ai meson exchanges as parts of the three- 
pion exchange. For r < 1 fm heavier meson exchanges 
and/or other degrees of freedom, such as quarks and glu- 
ons take place but the available theoretical calculations 
are not free from phenomenological parameters (see for 
instance Refs. 15] and 16]). Consequently, following the 
choice made in the case of the Paris NN potential [3] , 
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FIG. 2: Differential cross sections and polarization for the pp systems. The references to the experimental data can be found 
in Ref. 0. 



we use here an empirical short range real potential. 

As in Ref. [3| , we expand for r < 1 fm the phenomeno- 
logical radial potentials in power of r [Eqs. (|A3p and (|A4[) ] 
and match them to the theoretical ones at two point in 
the vicinity of 1 fm. Then, above 1 fm the theoreti- 
cal potentials are entirely preserved. For each isospin 
state, the spin structure of the NN interaction requires 
five independent invariants with five radial potentials. A 
phenomenological cubic expansion is used for the cen- 
tral components [see Eq. (|A3p ] and a quadratic one [see 
Eq. (jA4p ] for the other terms. This leads to nine param- 
eters representing the strength of the different empirical 
potentials at two (for the central components) or one 
(for the other) specific r values smaller than 1 fm (see 
Table |I| . The expressions of the full real potentials and 
more details on our fitting procedure are given in the 
Appendix. 

The imaginary potential Wj^fj includes the NN anni- 
hilation into mesons in the s-channel and can be calcu- 



lated from annihilation diagrams 0, [iBj ■ It implies the 
exchange of a baryon-antibaryon pair in the crossed t- 
channel and the resulting potential is then non-local and 
short range. We give in the Appendix the full expression 
we used. More details of its derivation can be found in 
the Appendix of Ref . 

The parameters of the real and imaginary short range 
parts are then determined by fitting the existing experi- 
mental data, viz. 915 data points in 1982, 3800 in 1994 
and over 4000 in 1999. References and discussions on the 
experimental 1999 data set is given in Ref H. A recent 
review on the NN data and the underlying physics can 



also be found in Ref [U 
used in the model 1999 



Our data consist of 4259 data 
plus the 64 data of the total 
np cross sections [9j and the ten level shifts and widths 
of the antiprotonic hydrogen of Refs. [l3| and [ll|. The 
best fit to this data set yields a x^/data = 4.52 for the 
present updated NN Paris potential. Here, the xVdata 
are calculated as in Ref. . The Paris 99 version [sj] has 
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FIG. 3: Total and differential cross sections and polarization for the pp — > nn systems. The references to the experimental data 
can be found in Ref. 0]. 



a corresponding x^/data = 4.59. The values of the 15 
parameters obtained in the present fit are compared to 
those of Paris 99 in Table HI 



2.2. Scattering observables 

In more detail, the value of x^/data are 6.66(9.76) 
for the 106 pp total cross section data (see Fig. [Ha)), 
2.66(4.47) for the 48 pp annihilation cross section data 
(see Fig.[IJb)), 2.03(4.29) for the 64 recent np total cross 
section data of Ref @ (see Fig. \V[c)), 0.87(1.56) for 
the 46 np annihilation cross section data (see Fig. [TJd)), 
2.29(3.08) for the 28 pp backward elastic differential cross 
section data of Alston Garnjost et al. (see Ref [6] in 
Ref and Fig.^a)), 4.17(4.05) for the 3392 pp elastic 
data (see more elastic data in Figs.[21[b-d) and 6.74(6.84) 
for the 639 pp —>■ nn charge exchange data (see some ex- 
amples in Fig. [3]). The x^/data are 6.71(11.75) for the 



ten level shifts and widths for the antiprotonic-hydrogen 
data[i3,liil (seeTablcini. All the above /data quoted 
in parenthesis correspond to calculations with the Paris 
99 potential [8j . Significant improvements were obtained 
with the Paris 99 potential compared to the earlier Paris 
82 and 94 [tI versions. As seen from the above de- 
tailed x^ and in comparison to the Paris 99 version 0], 
the present potential yields an improved description of 
all experimental observables but the pp elastic data. 

The results of the fit to the total and annihilation pp 
and np cross sections are shown in Figs, l(a-d). For 
Tiafc ^100 MeV the present model is closer to the data 
than the Paris 99 potential. The recent np total cross 
section (Fig. [Ijc)) is well reproduced. One has (see 
Fig. EUa)) a better description of the backward elastic 
differential cross sections for T^ab ^130 MeV. A sam- 
ple of pp elastic differential cross section at 288.3 MeV 
(Fig.H^b)) and polarization at 146.3 (Fig.H^c)) and 219.9 
MeV (Fig. [2Id)) are displayed in Fig. O As mentioned 
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TABLE II: Results of the fit to the level shifts Ai5_t and width of antiprotonic hydrogen for the present work compared 
to the experimental data. Units of AEl and Fl are keV for S waves and meV for P waves. Results for the Paris 99 
potential [3| are predictions. The corresponding Coulomb corrected scattering lengths are obtained from Eq. ([2]) for S-waves 
with a principal quantum number n = 1 and from Eq. for P-waves (n — 2). One has a{pp) = [a(T=0) + a(T=l)]/2, 
a(S'-world) = [a(singiet) + 3a(triplet)]/4 and a(Sum-P) = [3a(^Pi) + Sa^Pi) + 5a(^P2)]/ll. We use here the standard 
spectroscopic notation ^^^Lj for a given partial wave of spin S, of angular momentum L and total angular momentum J. 
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TABLE III: Binding energy in MeV of the close to threshold 
quasibound states of the present model and of the Paris 99 
potential 
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TABLE IV: Close to threshold resonances of the present 
model. The numbers in parenthesis correspond to the '^Pi 
resonance of Paris 99. The ^^Po and ^^Pi resonances have 
identical positions in the Paris 99 model. There is no ^^Pq 
resonance in the Paris 99 potential. 
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above, the present model is not as good as the Paris 99 
version on pp elastic data. Despite a slight improvement 
on the total x^/data for the pp —> nn charge-exchange, 
CEX, data (see above), it can be seen in Fig. [3] that both 
potentials fail to give a good fit. Had we allowed a 0.85 
normalization factor the reproduction of the integrated 
CEX cross section (see Fig. ^a)) would have been bet- 
ter. Note that at 147.3 MeV (Fig. [Sib)) a normalization 
factor of 0.74 is needed to reproduce the data. None of 
the two versions fit well the CEX polarization data as 
seen for instance in Fig. (ni^c) and (d). 



2.3. Antiprotonic hydrogen level shifts 

A comprehensive description on the protonium is dis- 
played in Ref. [T5|. Here, the S'-wave atomic level shifts 



AEs and widths Ts are related to the Coulomb corrected 
5-wave pp complex scattering length af by (see Eq. (3.7) 
of Rcf. 171), 
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.Fc 



27r 
f^pp 



1 - 3.154 



B 



(2) 



This relation is accurate to second order in af / B, good 
enough for the pp system. In Eq. ^ fipp = Mp/2 is the 
pp reduced mass, Mp — 938.27 MeV being the proton 
mass. The Bohr radius is B = l/{aiJ,pp) = 57.6399 fm, 
where a = 1/137.036 is the fine structure constant. 
The modulus of the zero range Coulomb wave function 
^coul(Q) gjyg^ 7/,coul(-o) ^ 1/{ttB^). It can be 

seen that for a scattering length of the order of 1 fni, 
the second term of Eq. ([2]) is a correction of a few per 
cent. Such a correction is negligible in higher angular 
momentum states and for P-wave (n = 2) we use [18| 



AEp - 



Tp 



IGfippB^ 



(3) 



We follow Eqs. (4) to (6) of Ref. [3 to calculate the 
Coulomb corrected scattering lengths. The phase shifts 
are obtained by solving the pp Schrodinger equation in 
configuration space with the Paris NN optical potential 
plus the Coulomb potential. 

The resulting S- and P-wave antiprotonic-hydrogen 
level shifts and widths are compared to the experimental 
ones [l3, [Hi in Table [Hi There is an overall improve- 
ment in comparison with the predictions of the Paris 99 
potential. Furthermore, our predictions of ^Pi, ^Pi and 
waves are given in this table. We also list all the 
corresponding Coulomb corrected scattering lengths re- 
lated here to the level shifts and widths through Eqs. ([2]) 
and ^ for S and P waves, respectively. Some of these 
numbers can be compared to those given in Table III of 
Ref. [iBl , where a NN potential is derived from a quark- 
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FIG. 4: Resulting central real ?7jvjv('", 0) and imaginary Wjvjv(^iO) potentials [see Eq. ([T])] compared with those of Paris 99 
The detailed definition of these potential can be found in the Appendix. 



based TV TV model Q supplemented by a state indepen- 
dent, phenomenological imaginary potential of Gaussian 
type. Energy shifts and widths show differences. All 
models give a value of AE{'^Po) about half of the exper- 
imental value. Due to the dominance of the one pion 
exchange for the ^Pi state, these models have similar re- 
sults for AE{^Pi) and T{^Pi). 



2.4. Bound states and resonances 

We search for the close to threshold bound states or 
resonances present in our model and in the Paris 99 po- 
tential Q. As seen in Figs. [4^ and[4|:, the real central sin- 
glet and triplet potentials have relatively strong medium 
and short range attractive parts. If the imaginary poten- 
tials are set to zero, several bound or resonant states ex- 



ist, however many of them disappear when the necessary 
annihilation is introduced [2l|. Following the method of 
this last reference, we find an isospin T=0 S-wave quasi- 
bound state "5*0 (using the notation ^'^+^ ^^~^^Lj) of 4.8 
MeV binding and of 52 MeV total width (see Table HIT)) . 
Such a state is absent in the Paris 99 [1| , Paris 94 Q and 
Paris 82 ^ potentials. A 9 MeV bound ^^Pi state of 18 
MeV total width is found. This triplet P-wave state is 
also found in the Paris 99 (see Table [111)) and Paris 94 (see 
Table VI of Ref. 0) models. The well known "p^ res- 
onance, resulting from an attractive one-pion exchange 
force, is present in all potentials (see Table [TV)) . This 
state has also been found in the recent NN constituent 
quark model of Ref. [3 ■ These close to threshold quasi- 
bound states and resonances are very difficult to detect 
in NN scattering data as, right from threshold, many 
partial waves contribute (see for instance the Paris 94 
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FIG. 5: As in Fig. |4] but for the spin-orbit and tensor optical potentials. 



phases in Ref. Chances to observe these states are 
larger in reaction processes which select the partial wave 
contributions such as in J/vl/ decays [sl. [22j. 



2.5. The optical potential 

The values of the short-range parameters of the poten- 
tial listed in TableUtogether with the plots in Figs.[4lja-d) 
of the central optical potentials at threshold (T^ab = 0) 
show the differences between the present potentials and 
those of 1999 [8]. The isospin real central singlet poten- 
tials are very similar (see Fig. |3Ia)). The isospin 1 real 
central singlet potential, less attractive for r > 0.5 fm 
than that of Paris 99, has a much deeper short range part 
(see Fig. IH^a)) and compare the corresponding C/q (''s) 
and Ug{r2) in Table H The imaginary T = central 
singlet potential is more repulsive than that of Paris 99 



(see Fig. SKb)) and its energy dependence is different as 
can be seen from the corresponding / values in Table [H 
The present imaginary T — I singlet potential, close to 
that of Paris 99 (see Fig.dfb)) is weak, however both po- 
tentials differ in their energy dependence (compare their 
/ values in Table [I]). The real central triplet potentials 
are quite attractive (see Fig. \^c)) which supports the 
•^P-wave quasibound states and resonances listed in Ta- 
bles ini and III. The imaginary T = and T ^ 1 triplet 
potentials are more repulsive than those of Paris 99 (see 
Fig. m^d)). The short range real spin-orbit potentials 
(Fig-EIa)) are very similar in both models while those of 
the tensor are stronger in the present model (Fig. [SKc)). 
The imaginary spin-orbit (Fig.[5Jb) and tensor (Fig.[5Kd) 
are much smaller (compare also the gLS and gT given in 
Table U). 
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3. DISCUSSION 

As can be seen in Figs. |3^, Ht, and[SJ: there is a more 
or less sharp change at the matching radius r = rc with 
Tj, = 0.84 or 1 fm (see the Appendix) between the phe- 
nomenological part of the potential and the theoretical 
part deduced from the Paris NN potential. This varia- 
tion is not due to the real part of the annihilation po- 
tential which in our model is very short ranged. As 
reminded in the Appendix, the NN annihilation into 
mesons, with intermediate nucleon-antinucleon state in 
the crossed channel, has been shown by dispersion tech- 
niques in Ref. [7] to lead to a very short ranged absorptive 
potential. The real part due to annihilation is expected 
to be also very short ranged. The high-partial wave anal- 
ysis of the NN interactions performed in Ref. [Is'l with 
the Paris model shows indeed that there is some limit of 
applicability for the G transformed potential. The uncer- 
tainty of the meson exchange interaction increases while 
its range becomes smaller. Then, in the 1 fm range and 
below, the real phenomenological short ranged part is at- 
tributed mainly to two sources: (i) to exchange of heavier 
mesons other than those explicitly included in the model, 
viz. TT, 27r, oj and Ai, (ii) to involved quark-antiquark 
(and gluon) forces based on QCD. The relatively sharp 
junctions between the real theoretical medium ranged po- 
tentials and the real phenomenological short ranged ones 
might then indicate some of these uncertainties in the 
G-parity transformed meson-exchange forces. 

The present updated Paris potential gives a better de- 
scription of all observables but the pp elastic data. As 
noted above in Sec. 12.21 it has only a slightly better 
overall x^/data of 4.52 versus 4.59 (Paris 99) for the 
4333 data considered. However the low energy data, 
TLab S 50 MeV, including the atomic data are better de- 
scribed (see below). We do not pretend that the present 
model is a better one but it is a different and an interest- 
ing one compared to the 1999 version in the sense that it 
generates a close to threshold quasibound ^^5*0 state. 

It has to be stressed that the actual contribution of 
the -^^Sq partial wave to the cross sections is small by the 
statistical reasons. In this sense the parameters of this 
wave cannot be determined precisely. This statement is 
mainly valid for TLab ^ 50 MeV when the P-wave con- 
tributions start to be very large. The atomic level shift 
and width in the spin singlet state are more constraining. 
However, even there, the contribution of the spin-0 state 
is only 1/3 of the contribution from the spin-1 state and 
the total contribution of the ^^Sq wave to the spin sin- 
glet atomic shift is about 25%. Fortunately, the new hp 
scattering data (which involves T = I amplitudes) allow 
one to fix the isospin 1 contribution in a better way. In 
particular the absorptive parts of scattering amplitudes 
are larger. It is the consistency of scattering data and of 
the atomic widths that induce the energy dependence of 
absorptive potential Wj^p;f{r, TLab) to differ from the pre- 
vious models. This reflects upon the position of the ^^5o 
state. In all Paris models, the potential in ^-^Sq state has 



TABLE V: The spin-singlet (not corrected for Coulomb) scat- 
tering lengths of the present model and of the Paris 99 poten- 
tial 





Present work 


Paris 99 




0.684-i 0.473 


0.979 -iO.294 




1.115-i 0.856 


0.844-i0.329 



been strongly attractive. Deeply bound states have been 
formed with binding energies in the few hundred MeV 
range. These states are not realistic as the extrapolation 
of the Wj^f[{r,TLab) so far down below the threshold is 
very uncertain. It is the next state which may or may 
not exists close to the threshold that is of our interest. 
It does not exist in Paris 99 version and arises in the 
present Paris one. The main difference that creates it, 
is the energy dependence of the Wj^]s}{r,TLab) generated 
by the hp data and the atomic level widths. 

For completeness, we show in Table FVl the spin-singlet 
(not corrected for Coulomb) scattering lengths in both 
models. The increase of absorption is visible, Im a(^^S'o) 
is determined fairly precisely as both Im a(^^S'o) and 
Im 0(^5*0) are now better known from the hp data and 
the antiprotonic hydrogen atom widths, respectively. It 
is not the case with the real parts which are larger than 
those obtained from atoms as seen in Table |TT1 On the 
other hand Figs. [T^ and [1)3 indicate that our total and 
annihilation pp cross sections are too small at Tiab ^ 50 
MeV. It shows some inconsistency, but it is also clear that 
the present potential is doing much better than Paris 99 
in this low energy and atomic region. This inconsistency 
does not necessarily reflect on the existence of the quasi- 
bound state. There are two reasons for that, one being, 
as mentioned above, the low statistical weight of the ^^Sq 
state. The other is the fact that Re a(^^iS'o) is not easily 
related to the bound state energy. Due to the peculiar 
shape of central singlet potential (Fig. IH^a) , the effective 
range expansion has a very short convergence radius and 
one cannot infer the binding energy from the scattering 
length. 

Let us also add that we have not been motivated to 
obtain the quasibound state. We do not sacrifice the low 
energy data at the expense of atomic data. The best 
fit gives an "in between" solution. Note that the partial 
values of for data limited to Tiab ^ 50 MeV with the 
atomic data (altogether 413 data) are 3544 for the present 
model and 5218 for the Paris 99 solution. A sizable frac- 
tion of these comes from charge exchange reactions. 
In fitting, we obtain several almost equivalent solutions 
as far as the minimum is concerned, although the 
present model has the best overall x^. These alternative 
solutions did not always produce the ^^Sq quasibound 
state. 
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4. SUMMARY AND CONCLUSIONS 

We have redetermined the short-range parameters of 
the Paris NN optical potential by fitting, besides the 
set of data used in the previous 1999 version, recent 
antiprotonic-hydrogen level shifts and widths and to- 
tal hp cross sections. Improvements in the fit are ob- 
tained. This model predicts quasibound states close to 
the threshold in the pp(^^So) and pp{'^^Pi) waves. Exis- 
tence of these states indicates a strong dependence on the 
parameters of the model. There is also a well established 
resonance in the pp(^^Po) wave. The ^Sq state repro- 
duces well [3] the recent jpp spectrum measured by the 
BES Collaboration [1]. The observed peak in the invari- 
ant mass of the produced mesons in the J/ip — + ^TT^n^rj' 
decay (2^ can be explained by an interference of the qua- 
sibound state ^^5*0 with a background amplitude [22 1 . 
Let us mention that for uses of the Paris NN potential, 
like for the initial or final state interactions in various 
processes, the present potential can be provided upon 
request. 
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APPENDIX A: FULL EXPRESSION OF THE 
OPTICAL POTENTIAL 

For completeness, the full expression of the optical 
potential is revisited below. For each isospin value 
T = or T = 1, the real potential can be expressed in 
terms of the five usual nonrelativistic invariants, flo = 
(l-*Ti •cr2)/4, rii = (3 cTi • ^2) /4, Qls = L-S, 
fir = 3((Ti ■ r (T2 ■ r)/r^ - cri ■ 0-2 and ^302 = [cri ■ 
L 0-2 • L + 0-2 • L (Ti • L)/2. It reads 

f/jVJV (r, TLab)) = Uo (r, TLab) ^0 + Ui (r, TLab) ^1 

+ ULs{r)nLs + UT{r)nT 

+ Uso2{r)^so2- (Al) 

The linear nonlocality in the central singlet, Uq and cen- 
tral triplet, Ui potentials is expressed as 



U{r,TLab)^U%r)+TLab U\r). 



(A2) 



The potential, for r > (r^ < 1 fm), is the G-parity 
transform of the theoretical Paris NN potential |13l] . We 



use, for the w exchange (m^^ = 782.7 MeV), g^/^n = 
11.75 as in Ref. [l^. However, in order to have a more 
attractive isospin central singlet potential in the vicin- 
ity of 1 fm, we modify the coupling of the shorter range 
Al exchange (m^^ = 1100 MeV) from = 14. to 10.4. 

For r < Vc, the empirical potentials are given by a cubic 
r expression for the central UqIt) and C/f (r) terms: 



U (r) = a^r^ + 02^^ + air + oq 



(A3) 



and by a quadratic one for, UQ i{r), ULs{r), Urir) and 



U{r) = 62?-^ + hir + bo. 



(A4) 



The parameters Ui (i = to 3) and bi {i = to 2) are 
determined (i) by matching to the theoretical potential 
a,t r ~ tq ^ Tc and r = ri = tq + Ar with Ar = 0.15 fm, 
(ii) choosing a phenomcnological height at r2 — 0.587 fm 
and at r = 0.188 fm (here only for U§i{r)). For all 
isospin-0 potentials Tc = 1 fm and for all isospin- 1 terms 
Tc = 0.84 fm except for U§{r) where = 1 fm. While 
solving the Schodinger equation we have regularized the 
tensor potential Uxir) at small r by multiplying it by 



Fir) 



{prf 



1 + [prf 



(A5) 



with p = 10 fm""'^. 

As shown in the Appendix of Ref. the imag- 
inary potential VF^^(r, Ttaj,), arising from nucleon- 
antinucleon annihilation into mesons, can be approxi- 
mated by a short range radial function proportional to 
effective phenomcnological couplings with a linear energy 
dependence for the central and spin-spin components. 
One writes 



9C (1 + fcTLab) 

gSS (1 + fssThab) (Ti ■ (T2+ gx^T 

Koi2mr) 



Arn^ r dr 



(A6) 



where the modified Bessel function Ko{2mr) is the 
Fourier transform of a dispersion type integral resulting 
from the calculation of the NN annihilation box diagram 
into two mesons with a nucleon-antinucleon intermediate 
state in the crossed i-channel One has 



ifo(2mr) 



dt'- 



-Vt'r 



(A7) 



In Eqs. (|A6P and (jA7p . m is taken to be quite close to 
the nucleon mass, m = 940 MeV. To avoid the singular 
behavior at r = 0, we regularize the central and spin-spin 
potential of Eq. (jA6| by multiplying them with 



G{r) = (1 - e-2™'')4 



(A8) 
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The imaginary spin-orbit and tensor potentials are mul- 
tiplied by 

H{r) = (1 - e-2™y. (A9) 

The values of the empirical real potentials at r = r2 
and ra and of the parameters gt {i = C, SS, T and LS) 
and fi (i = C, SS), determined by the fit, are displayed 
in Table HI These 15 parameters play a different role in 



the fit. The important parameters are 

(i) the six values of the central singlet and triplet, of the 
tensor and of the L ■ S components at r = r2 = 0.587 fm 
for the real part, 

(ii) the four couplings of the singlet and triplet central 
terms for the imaginary part. 

A fine tuning of the fit is then obtained by adjusting the 
five remaining parameters. 
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